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ABSTRACT 

I present a systematic study of gamma-ray flares in blazars. For this purpose, I propose 
a very simple and practical definition of a flare as a period of time, associated with 
a given flux peak, during which the flux is above half of the peak flux. I select a 
sample of 40 brightest gamma-ray flares observed by Fermi /L AT during the first 4 
years of its mission. The sample is dominated by 4 blazars: 3C 454.3, PKS 1510-089, 
PKS 1222+216 and 3C 273. For each flare, I calculate a light curve and variations of 
the photon index. For the whole sample, I study the distributions of the peak flux, 
peak luminosity, duration, time asymmetry, average photon index and photon index 
scatter. I find that: 1) flares produced by 3C 454.3 are longer and have more complex 
light curves than those produced by other blazars; 2) flares shorter than 1.5 days in 
the source frame tend to be time-asymmetric with the flux peak preceding the flare 
midpoint. These differences can be largely attributed to a smaller viewing angle of 
3C 454.3 as compared to other blazars. Intrinsically, the gamma-ray emitting regions 
in blazar jets may be structured and consist of several domains. I find no regularity 
in the spectral gamma-ray variations of flaring blazars. 
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1 INTRODUCTION 

Blazars constitute the most numerous class of permanent 
gamma-ray sources in the sky. In the 2-year Fermi Large 
Area Telescope (LAT) Source Catalog (2FGL; Nolan et al. 
2012), they constitute 67% of associated point sources. They 
are characterised by strong and chaotic variability, occa- 
sionally producing spectacular flares. These events have at- 
tracted a considerable attention, triggering many extensive 
multiwavelength campaigns that attempted to characterise 
the overall behaviour of individual objects {e.g., Vercellone 
et al. 2010; Abdo et al. 2010b; Agudo et al. 2011; Acker- 
mann et al. 2012; Hayashida et al. 2012). Unfortunately, the 
results of these campaigns vary widely, revealing relatively 
few consistent patterns that can be used to constrain the 
underlying physics of relativistic jets. 

The Fermi Space Telescope offers a unique capability 
in the observational astronomy. Its main detector, LAT (At- 
wood et al. 2009) , has a very wide field of view, of ~ 50° . For 
most of its scientific mission starting in August 2008, Fermi 
operated in the scanning mode, allowing a uniform cover- 
age of the whole sky every several hours. Therefore, almost 
every blazar flare that happened during the Fermi mission 
was observed from the beginning to the end with basically 
the same sensitivity. It is then straightforward to produce a 
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fairly complete and unbiased sample of gamma-ray flares of 
blazars, and to compare their properties within individual 
sources and between them. Such an approach has a poten- 
tial to clearly demonstrate the similarities and differences 
among individual blazar flares. 

In this work, I study a sample of 40 bright gamma-ray 
flares of blazars detected by Fermi /LAT during the first 4 
years of its mission. In Section 2, I propose a useful def- 
inition of a gamma-ray flare, and describe the process of 
selecting the sample. In Section 3, I present light curves of 
these flares, and study the distributions of their duration, 
time asymmetry, average photon index and photon index 
scatter. In Section 4, I discuss the general properties of the 
flares and their interpretation. I conclude in Section 5. 



2 SAMPLE SELECTION 

Any discussion of blazar flares immediately leads to a fun- 
damental difficulty - how to define a flare. The nature of 
blazars variability is stochastic (Abdo et al. 2010c), charac- 
terised by a very high duty cycle, and there is no consistent 
'background' flux level. It is thus difficult to draw a division 
between 'flares' and 'quiescence'. In fact, many researchers 
question the validity of using the term 'flare' in the context 
of blazar light curves. On the other hand, this term is intu- 
itive, as blazars spend very little time in their highest flux 
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Table 1. List of the brightest gamma-ray flares of blazars. MJDp Ca k is the moment of flux peak, -F pea k is the peak flux in units 
of 10 — 6 phs - 1 cm - 2 , ti is the flux doubling time scale, ti is the flux halving time scale, T = (t\ + *2) is the flare duration, k = 
(ti — t2)/(t\ +t2) is the time asymmetry parameter, (T) is the average photon index, and rms(T) is the photon index scatter. All times 
are in units of days. References: a - Abdo et al. (2011), b - Ackermann ct al. (2010), c - Tanaka ct al. (2011), d - Abdo et al. (2010b), 
e - Abdo et al. (2010a). 
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states (Tavecchio et al. 2010). My opinion is that the term 
is certainly useful, but for a systematic study a rigorous def- 
inition is necessary. 

Here, I propose a very simple definition, in which a flare 
is a contiguous period of time, associated with a given flux 
peak, during which the flux exceeds half of the peak value, 
and this lower limit is attained exactly twice - at the be- 
ginning and at the end of the flare. One can show that this 
definition does not allow any two flares to overlap. Every two 
flares must be separated by a flux minimum which is lower 
than half of the lower peak. A flare may include minor peaks, 
but they are not considered to be peaks of separate subflares. 
This definition is natural in a stochastic light curve with a 
power density spectrum (PDS) in the form of a power-law, 
as it is based on equal flux ranges on a logarithmic scale. It 



is practical, as it is very easy to determine local flux peaks 
and the associated lower flux limits. Also, it automatically 
provides two standard variability time scales often discussed 
in the blazar studies - the flux-doubling time scale and the 
flux-halving time scale. Figure 1 illustrates this definition in 
various real light curves. 

The process of selecting a flux-limited sample of blazar 
flares is divided into two steps. In the first step, a preliminary 
list of blazar active periods is extracted from the Fermi /LAT 
Monitored Source List 1 . It is a database of daily and weekly 
flux estimates for all sources exceeding the 0.1 — 300 GeV 
photon flux of 10 -6 phs _1 cm~ 2 . I select all daily photon 
fluxes during the first 4 years of the Fermi mission (MJD 

1 http : //f ermi . gsfc .nasa.gov/ssc/data/access/lat/msl_lc/ 
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Figure 1. Light curves of individual gamma-ray flares calculated from the Fermi /LAT data, presented in the order of descending peak 
photon flux density. For each flare, the upper panel shows the normalised photon flux density integrated in the energy range 0.1 — 300 GeV 
in 0.5-day time intervals with a 0.1-day step, and the lower panel shows the corresponding photon index. The horizontal axes show time 
in days measured relative to the flare peak. Data points belonging to the flare period are highlighted with colours. 



54682 - 56143) above 3 x 1(T 6 phs^crrr 2 for sources lo- 
cated away from the Galactic plane (|6| > 10°). For each 
blazar, these daily entries are grouped into active periods, 
which are separated by gaps of at least 10 days. If an ac- 
tive period contains two flux peaks separated by at least 
10 days and a flux minimum lower than the half of the 
lower peak flux or at least one-day gap, it is split into two 
active periods along the flux minimum. The active peri- 
ods for all blazars are sorted according to the decreasing 
peak flux, and the first 30 active periods are selected for 
the second step. This choice corresponds to a minimum flux 
peak of 4.5 x 10~ 6 s _1 cm~ 2 , and the preliminary list of ac- 
tive periods consists of 6 blazars: 3C 454.3, PKS 1510-089, 
PKS 1222+216, 3C 273, PKS 0402-362 and PKS 1329-049. 

In the second step, I perform the maximum likelihood 
analysis of the Fermi /LAT data. I use the standard analysis 
software package Science Tools v9r27pl, with the instru- 
ment response function P7S0URCE_V6, the Galactic diffuse 
emission model gal_2yearp7v6_v0, and the isotropic back- 



ground model iso_p7v6source. Events of the SOURCE class 
were extracted from regions of interest of 10° radius centred 
on the position of each blazar. The background models in- 
cluded all sources from the 2FGL within 15°, their spectral 
models are power-laws with the photon index fixed to the 
catalogue values. I checked the raw count maps for any new 
source not included in the 2FGL. In the case of PKS 1510- 
089, the source model included TXS 1530-131 (Gasparrini 
& Cutini 2011). I selected events in the energy range be- 
tween 100 MeV and 300 GeV. The spectra of all studied 
blazars were modelled with power-laws. I used overlapping 
time bins, with the shift between consecutive bins equal to 
1/5 of the bin lengths. This allows to avoid the dependency 
of the results on the 'phase' of the time bins and provides 
a natural interpolation of the light curve. While the consec- 
utive flux measurements are not entirely independent, this 
is not a problem for the subsequent analysis of the flare pa- 
rameters. 

Light curves calculated with different time bins reveal 
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Figure 1. (continued) 



different amounts of detail, which affects the classification of 
certain events as flares, and the flares parameters. Shorter 
time bins reveal more details and potentially shorter flares, 
but limited photon statistics determines the size of the sam- 
ple of flares that can be studied reliably. The very brightest 
gamma-ray flares (> 1(T 5 phs" 1 cm" 2 ) can be probed once 
per orbital period of the Fermi satellite (~ 1.5 h; Foschini 
et al. 2011). Here, I decide to use the time bin length of 
0.5 d, with the shifts of 0.1 d, always starting from a full 
modified julian day (MJD). 

Systematic analysis of all candidate active periods re- 
vealed that some of them consist of a few flares, so that 
the list of all flares revealed by this analysis exceeds 50 
entries. From this, I select the final list of the 40 bright- 
est flares, which corresponds to the minimum peak flux of 
7.1 x 10 -6 phs _1 cm -2 . This limiting peak flux value cor- 
responds roughly to the peak of a histogram of the loga- 
rithmic peak flux distribution of all identified flares. The 
final list of brightest gamma-ray blazar flares is reported in 
Table 1. For each flare, I calculate the flux doubling time 
ti, the flux halving time t2, and the total flare duration 



T = ti + ti. I also calculate the time asymmetry parameter 
k = (ti — t2)/(t\ + t2), which can take values between — 1 
and 1, with k — corresponding to a time-symmetric flare 
(ti = ti); the average photon index T (Ne tx -E _r ), and the 
root mean square (rms) of deviations of the photon index 
from the mean value. Individual flares will be referred to 
using their row number in Table 1. 

Very short flares (e.g. flare #3 — PKS 1510-089 at 
MJD 55854) indicate that the measured flux can change 
significantly over a single time bin shift of St — 0.1 d. 
This provides a practical estimate of accuracy in determin- 
ing the moments of flare beginning, peak and ending. The 
uncertainty of parameters t\ , t2 and T can be estimated as 
\/~28t = 0.14 d. The uncertainty of the asymmetry parame- 
ter k is given by 5k = [2(3 + k 2 )] 1/2 {St/T), which is large for 
the shortest flares. The values of Sk are reported in Table 1. 



3 RESULTS 

My sample of the brightest gamma-ray flares of blazars re- 
ported in Table 1 is produced by only 5 blazars: 3C 454.3 
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Figure 2. Left panel: distribution of the observed flare duration T vs. the peak flux density -F pGa k- Right panel: distribution of the 
source-frame flare duration T/(l + z) vs. the peak apparent luminosity L pea it. The shape and colour of each point indicates the host 
blazar. 
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Figure 3. Distribution of the source-frame flare duration T/(l + 
z) vs. the time asymmetry parameter k. The shape and colour of 
each point indicates the host blazar, and the point size is propor- 
tional to the logarithm of peak flux -F pca k- 



(16 entries), PKS 1510-089 (14 entries), PKS 1222+216 
(6 entries), 3C 273 (3 entries), and PKS 0402-362 (1 en- 
try). Remarkably, 3C 454.3 produced 7 of the 10 bright- 
est flares. About 20 flares exceed the peak flux level of 
10~ 5 phs^cm- 2 . 

In the following analysis, I will compare both the ob- 
served parameters, and parameters measured in the source 
frame. The redshifts and luminosity distances of the 5 
blazars that produced the sample are given in Table 2. 



Table 2. Redshifts and luminosity distances of blazars produc- 
ing the sample of brightest gamma-ray flares. I adopt standard 
ACDM cosmology with H = 71 kms -1 Mpc -1 , H m = 0.27 and 
fl A = 0.73. 



blazar 


2 


rJ L [Gpc] 


3C 454.3 


0.859 


5.49 


PKS 1510-089 


0.360 


1.92 


PKS 1222+216 


0.432 


2.37 


3C 273 


0.158 


0.75 


PKS 0402-362 


1.423 


10.31 



3.1 Light curves 

Figure 1 shows normalised light curves and photon index 
variations of all 40 flares, in the order of decreasing peak 
flux, the same as in Table 1. The data points belonging to 
the flare period are highlighted with a colour. One can im- 
mediately notice a great variety of the flare shapes and du- 
rations. Looking only at the first 10 panels, it is striking 
that the flares produced by 3C 454.3 appear very different 
from those produced by PKS 1510-089 and PKS 1222+216 

— the former are longer and more complex than the latter. 
However, looking at the next 30 panels, the picture is less 
clear. Many flares, especially the longest ones, are charac- 
terised by a complex structure with multiple peaks separated 
by shallow minima. The time scale of these substructure is 
comparable to the time scale of short flares characterise by 
a single peak. Medium-long flares (T ~ 2 d) often show a 
single substructure following the main peak (e.g., flare #12 

- PKS 1222+216 at MJD 55366), and sometimes preceding 
the main peak (e.g., flare #20 — PKS 1510-089 at MJD 
54917). This results in the significant time asymmetry of 
these flares. 



© 0000 RAS, MNRAS 000, 000-000 



6 Nalewajko 




1 0.5 



1 0.5 



1 0.5 



Figure 4. Scatter plots showing the normalised 0.1 — 300 GcV photon flux (horizontal axes) vs. the corresponding photon index (vertical 
axes) for individual flares. The intensity of data points indicates the time flow from dark to bright. 



3.1.1 Flare duration 

In order to further investigate the temporal properties of 
blazar flares, in Figure 2, I plot the distribution of the flare 
duration vs. peak flux. Flare durations span the range be- 
tween 0.5 and 10 days. The lower limit corresponds to the 
time bin length, and the number of flares with T < 1 d (9) 
is the same as that of flares with 1 d < T < 2 d. While the 
sample might be somewhat incomplete for T < 1 d, it is not 
likely for longer flares. The average flare duration is 2.7 d, 
and the median is 2.1 d. There appears to be no general 
trend between the flare duration and the peak flux. How- 
ever, there is a significant difference between 3C 454.3 on 
one hand, and PKS 1510-089, PKS 1222+216 and 3C 273 
on the other. Flares produced by 3C 454.3 are on average 
longer than those produced by other blazars. 3C 454.3 did 
not produce a flare shorter than 1.6 d, while other blazars 
generally produced flares shorter than 3 d, with the sole ex- 
ception of flare #28 (PKS 1510-089 at MJD 55991). For 
the peak flux greater than 1.2 x 10 -5 phs _1 cm -2 , there 
is a complete separation of the two groups. Also, there is 
no discernible difference between the distributions of flares 
produced by PKS 1510-089, PKS 1222+216 and 3C 273. 
Also in Figure 2, I plot the distribution of the source- 



frame flare duration vs. peak apparent luminosity. I do not 
use luminosity integrated from the power-law spectral fit 
in the 0.1 — 300 GeV range, because: 1) the correspond- 
ing source-frame energy range depends on the source red- 
shift, 2) the gamma-ray spectra of these blazars are actually 
curved /broken, so the power-law spectral fits overestimate 
the bolometric luminosity above a few GeV. Instead, I use 
the vL v value corresponding to the source-frame photon en- 
ergy of 500 MeV. Apparent luminosities of the brightest 
blazar flares span three orders of magnitude, mainly due to 
the differences in luminosity distances between the sources. 
The flare produced by PKS 0402-362, while 39th in observed 
photon flux, is 3rd in apparent luminosity, due to the high 
source redshift. There is no clear trend between flare lumi- 
nosity and duration. 

3.1.2 Flare asymmetry 

Further insight into the flare profiles can be gained by look- 
ing at the time asymmetry parameter k. In Figure 3, I plot 
the distribution of the source- frame flare duration T / (1 + z) 
vs. k. The values of k range between -0.7 and 0.7. This means 
that the observed flares do not show a very rapid onset di- 



© 0000 RAS, MNRAS 000, 000-000 



The brightest blazar flares 7 




0.5 1 0.5 1 0.5 1 0.5 1 0.5 1 



Figure 4. (continued) 



rectly to the peak or a very rapid decay directly from the 
peak. The average value of k is — O.f, and the median is 
—0.16. Therefore, there is a slight tendency for the flux dou- 
bling time ti to be shorter than the flux halving time ti. 
Figure 3 reveals a significant difference between short and 
long flares. Flares shorter than ~ 1.5 days in the source- 
frame have predominantly negative k, while flares longer 
than ~ 1.5 days have a broad and roughly uniform distri- 
bution of k. Individual blazars produce flares of different 
distributions of k, which can be largely understood as re- 
sulting from different distributions of the flare durations. 
3C 454.3 produces flares of a broad distribution of k, with 
4 out of the 5 shortest flares characterised by negative k. 
PKS 1510—089 produces flares with mostly negative k, how- 
ever, the 3 longest flares have significantly positive values of 
k. PKS 1222+216 produces fairly symmetric flares, while the 
3 flares produced by 3C 273 have significantly negative k. 
The most asymmetric flare in the sample is the one produced 
by PKS 0402-362. There is a slight trend for the brighter 
flares to have more negative values of k; the average value 
of k for the 20 brightest flares is —0.16. 



3.2 Spectral variations 



In order to study variations of the photon index in detail, 
in Figure 4, I show scatter plots of the normalised flux 
density F pe ak vs. the photon index F for each flare in the 
sample. Varying point intensity is used to show the flow of 
time. It is clear that individual flares are characterised by 
different average values of the photon index, and also by 
different spreads of its values. Scatter plots for long flares 
are obviously more complex than those for short flares. In 
some cases, the spectral evolution of the gamma-ray flare 
takes form of a clear loop (e.g. flare #11 — PKS 1510-089 
at MJD 55868). Such loops were observed numerously in 
blazars in gamma rays (e.g., Ackermann et al. 2010) and X- 
rays (e.g., Fossati et al. 2000). The most basic characteris- 
tics of a loop is its orientation in time. In Figure 4, clockwise 
(CW) loops correspond to spectral hardening (decreasing V) 
at the peak, and counterclockwise (CCW) loops correspond 
to spectral softening (increasing F) . I can identify clear cases 
of CW loops (flares #8, #9, #11) and CCW loops (flares 
#23, #26, #30). However, most cases are difficult to clas- 
sify, and thus I cannot draw a definite conclusion about the 
statistics of loop orientations. 
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3.2.1 Average photon index 

In Figure 5, I plot the distribution of the source- frame flare 
duration T/(l+z) vs. average photon index (F). The average 
photon index takes values between 1.9 and 2.6. The extreme 
values are realised only by very short flares, with (1 + z)T < 
1 d, and for long flares the range of (F) converges to the 
vicinity of 2.3 - 2.4. The average value of (F) for the sample 
is 2.27, and the median is 2.3. Individual blazars produce 
flares of significantly different distributions of (F) . 3C 454.3 
produces flares characterised by 2.2 < (r) < 2.45, i.e. with a 
hardness typical for the whole sample and with a very small 
scatter. PKS 1510—089 produces flares with a wide range 
of average photon index, 2 < (r) < 2.6. PKS 1222+216 
produces relatively hard flares, with 1.9 < (F) < 2.25, and 
3C 273 produces relatively soft flares, with (F) ~ 2.4. There 
is no clear dependence between the peak flux density Fpeak 
and the average photon index. 



3.2.2 Photon index scatter 

In Figure 6, I plot the distribution of the average photon 
index (r) vs. the photon index scatter rms(F). The maxi- 
mum value of the photon index scatter is ~ 0.13, and the 
minimum value is ~ 0.02. The average and median values 
are ~ 0.09. I find the general trend that the photon in- 
dex scatter is higher for softer flares (higher (r)), the Pear- 
son's correlation coefficient is p((r) ,rms(r)) = 0.5. Individ- 
ual blazars show a large range of rms(r) values. Flares of 
PKS 1222+216 are characterised by low photon index scat- 
ter (rms(F) < 0.1), while flares of 3C 273 show high average 
(rms(r) > 0.09). Flares of duration longer than 0.9 d in the 
source frame generally have rms(F) > 0.05, but the high 
photon index scatter values (rms(r) > 0.1) are exhibited by 
both long and short flares. 



4 DISCUSSION 

Variability of blazars is a stochastic phenomenon, regardless 
of in which band and at what time scales it is studied. Even 
the very notion of a flare, while intuitive, is controversial, 
because it is difficult to formulate a simple and practical 
definition. Here, I propose a simple mathematical definition 
of a flare, based on the peak and half-peak flux levels, with- 
out considering the complexity of the observed light curves. 
Such defined flare may be fairly simple, with a single roughly 
time symmetric peak, or complex, with several local maxima 
and minima. In practice, even such a simple definition of a 
flare faces some problems: the flare duration and peak flux 
depend on the sampling (in the case of Fermi /LAT data, 
the time binning), and limited photon statistics introduces 
some ambiguity in determining the flare beginning and end. 
Nevertheless, using the same definition for a uniform sample 
of blazar flares reveals new information about the variability 
of blazars. 

The focus of this study is on individual flares hosted 
by a handful of blazars. Flares of any given blazar, sepa- 
rated by months or years, can be considered independent 
phenomena, as they are produced in constantly flowing rel- 
ativistic jets. Any similarities and differences between them 
reflect independent realisations of the underlying processes 
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of energy dissipation and particle acceleration. Similarities 
and differences between flares produced in different blazars 
reflect the environmental factors that vary over time scales 
of many years. 



4.1 Temporal properties 

The most appealing finding of this study is the systematic 
difference between flares produced by 3C 454.3 and other 
bright blazars. Flares produced by 3C 454.3 are significantly 
longer and more complex. A brief look at the first 10 pan- 



© 0000 RAS, MNRAS 000, 000-000 



The brightest blazar flares 9 



els of Figure 1 shows that they are characterised by mul- 
tiple secondary maxima and minima. In accordance with 
the adopted definition of the flare, these minima are always 
higher than half of the peak flux for a given flare. The light 
curve of 3C 454.3 can be contained in a limited range of flux 
values for several days. Other blazars can produce such long 
flares only occasionally (e.g., flare #28 in PKS 1510-089). 
The long flares of 3C 454.3 can be described as a collec- 
tion of closely-spaced simple flares of comparable peak flux 
levels. On the other hand, short flares displayed mostly by 
other blazars appear to be simple (3C 454.3 can also produce 
simple flares, e.g. flare #31). 

The distinct structure of flares produced systematically 
by 3C 454.3 could be explained by its special jet orientation. 
Hovatta et al. (2009) provide estimates of the jet Lorentz 
factor Tj = (1 — /3 2 )" 1 / 2 , where /3j = Uj/c, and viewing angle 
#obs for a large sample of blazars monitored in the MO JAVE 
program. They estimate the viewing angle of 3C 454.3 at 
Sobs — 1-3° — 0.45r~ , which is significantly smaller than 
the values for PKS 1510-089 (3.4°), PKS 1222+216 (5.1°) 
and 3C 273 (3.3°). The corresponding Doppler factor for 
3C 454.3, V = [rj(l - Pjcose^)}- 1 ~ 33 ~ l.7T u is the 
highest among the blazars producing the brightest gamma- 
ray flares. This explains why 3C 454.3, at z — 0.86, being 
more distant than other blazars, is the brightest of them. In 
fact, the viewing angle of the inner jet of 3C 454.3 could be 
even smaller than that estimated from the VLBI studies, as 
these estimates by their nature are not very precise, and the 
jets may be slightly bent even at parsec scales (Hough et al. 
2002). The observed emission from a jet element propagating 
at the Doppler factor T> is boosted by factor T> 4 . For o bs = 
we have T> = ( 1 + /3j ) T j , while for # ODS = rr -1 (Doppler cone) 
we have T> = Tj. The difference in the Doppler factors by 
factor (l+/3j) ~ 2 corresponds to a difference in the observed 
fluxes by factor ~ 16. For an average blazar, one can expect 
that # bs ~ 0} ~ ly 1 , where 9j is the jet half-opening angle, 
so that the observer is located at random orientation within 
the Doppler cone 2 . In such a case, only a fraction of the 
jet volume is strongly Doppler-boosted and effectively con- 
tributes to the observed emission (see Figure 7). However, 
when o b s <C dj, most of the jet volume is strongly Doppler- 
boosted. If individual simple flares are produced in different 
parts of the jet cross-section, their Doppler factors will be 
very similar, resulting in similar peak fluxes. Thus, the long 
complex gamma-ray flares of 3C 454.3 are consistent with 
its very close orientation to the line of sight. 

Another interesting finding is that short flares (T < 2 d) 
tend to be time-asymmetric, with the flare onset shorter 
than the flare decay. In some cases, it can be seen that 
strongly asymmetric flares consist of a main peak and a 
subsequent secondary maximum. In the case of flare #12 
in PKS 1222+216, a hierarchical structure of simple flares 
is present. Why should secondary maxima follow the main 
peak, rather than preceding them? It can be understood 

2 Pushkarcv ct al. (2009) estimated the jet opening angles for the 
MOJAVE blazars. Their results indicate that FjOj < 1, which 
would not allow for high contrasts in the Doppler factors. How- 
ever, this measurement, made at the frequency of 15 GHz, probes 
jet regions located significantly beyond the gamma-ray emission 
sites, as indicated by a systematic delay of the 15 GHz radio 
signals with respect to the gamma rays (Pushkarev et al. 2010). 




Figure 7. Effects of the viewing angle o bs on the relative 
Doppler boost (darker colour - higher Doppler factor T>) of local 
emitting regions in a conically expanding relativistic jet. Observer 
1 is located at o bs "C 0j, and observes similar Doppler factors 
from emitting regions located at different azimuthal angles. Ob- 
server 2 is located at o b s ~ 9j (~ Tj" 1 ) ant l observes different 
Doppler factors, and radiation from a more strongly boosted emit- 
ting region is observed earlier because of shorter light travel time. 
I propose that Observer 1 corresponds to us watching 3C 454.3, 
and Observer 2 corresponds to us watching other blazars produc- 
ing the sample of bright flares. 



in the scenario discussed in the previous paragraph. Emit- 
ting regions located at the same distance measured along 
the jet, but oriented at smaller angles to the line of sight, 
are located closer to the observer, so the light travel time is 
shorter. At the same time, they are more strongly Doppler 
boosted, therefore one observes a major peak followed by 
minor peaks. In the case of 3C 454.3, where # b s <C dj, the 
Doppler factors of different emitting regions are similar, and 
thus its flares appear to be time-symmetric. 

If the different jet orientations between 3C 454.3 and 
other blazars explain the major differences between param- 
eters of flares produced by these sources, I can expect that 
blazar flares are intrinsically similar in all sources. The sce- 
nario presented in Figure 7 assumes that the gamma-ray 
emitting regions in blazar jets are structured and consist of 
several domains. The long multiply-peaked flares of 3C 454.3 
indicate that many simple flares of similar peak fluxes can 
be produced within a short range of time, so that they par- 
tially blend with each other. This partial blending may indi- 
cate that these simple flares are not entirely independent. A 
possible scenario for their origin is a large-scale jet pertur- 
bation that is intrinsically unstable or propagates through 
an inhomogeneous medium. Many major blazar flares have 
been associated with disturbances propagating along the jet 
that later-on are identified with superluminal radio knots 
(Marscher et al. 2012). Such disturbances can be interpreted 
as shock waves if the jet is weakly magnetised (Lind & 
Blandford 1985), or reconnecting magnetic domain inter- 
faces if the jet is highly magnetised (Lovelace et al. 1997; 
Nalewajko et al. 2011). As most theoretical models for the 
formation of relativistic jets predict that they are initially 
significantly magnetised (Blandford & Znajek 1977; Bland- 
ford & Payne 1982), and that it is difficult to efficiently 
convert the magnetic energy flux into kinetic energy flux 
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(Begelman & Li 1994), scenarios based on magnetic recon- 
nection appear to be more likely. Even if the magnetic do- 
mains cannot be easily reversed in AGN jets, reconnection 
could be triggered by short-wavelength current-driven insta- 
bilities (CDI; Begelman 1998; Nalewajko & Begelman 2012). 
On the other hand, shock waves can be efficient in magne- 
tised jets if the jet features highly evacuated gaps (Komis- 
sarov 2012). 

Variability of bright Fermi /LAT blazars during the first 
year of the Fermi mission was studied comprehensively by 
Abdo et al. (2010c). In particular, they studied individual 
flares of 10 sources of high variability index. They used 
light curves integrated using 3-day bins and fitted them 
with double-exponential flare templates allowing for arbi- 
trary time asymmetry. They found that the average flare 
duration was (T) ~ 12 d and the average flare asymmetry 
(k) ~ 0.08, with broad distributions of these parameters. 
While both parameters were calculated in a different way 
than in this work, the flares studied by Abdo et al. (2010c) 
are different events than the flares studied here. By using 
time bins of 3 days and looking at entire light curves, they 
picked many long gently varying components. In this work, 
I use 0.5-day time bins and look only at pieces of the light 
curves that belong to individual flares according to my defi- 
nition, and I find no flares longer than 10 days. While Abdo 
et al. (2010c) conclude that no particular type of flare asym- 
metry (positive/negative) is being favoured by their results, 
I show that only flares shorter than ~ 1.5 d tend to be 
asymmetric. 

4.2 Spectral properties 

Average gamma-ray spectral index of the brightest blazar 
flares ranges between 1.9 and 2.6. As all of the studied 
blazars belong to the subclass of flat-spectrum radio quasars 
(FSRQs), these values are rather typical (Ackermann et al. 
2011). I note that for the longest flares spectral indices con- 
verge into a narrow range between 2.3 and 2.4. This is most 
likely a statistical effect, as the average is calculated from 
a larger number of data points. Long flares show complex, 
chaotic structures in the gamma-ray flux - photon index 
space, and the scatter of photon index values is of order 
±0.1. But even short flares, while they usually have simple 
light curves, often show a substantial photon index scatter. 
Individual blazars show significant differences in the spectral 
behaviour. PKS 1222+216 typically shows relatively hard 
spectra with small photon index scatter, while 3C 273 shows 
relatively soft spectra with high photon index scatter. Over- 
all, the gamma-ray spectral variability of flaring blazars is 
completely irregular. Because of lack of consistent patterns 
(e.g. CW or CCW loops), the flux - photon index diagrams 
are not very useful as probes of the physics of particle ac- 
celeration in FSRQs. 



5 CONCLUSIONS 

I selected a sample of the 40 brightest gamma-ray flares of 
blazars that occurred during the first 4 years of the Fermi 
mission. A flare is defined as a period of time over which the 
flux exceeds half of the peak value. For each flare, I calculate 
its peak flux, peak luminosity, duration, time asymmetry, 



average photon index and the photon index scatter. I present 
light curves and flux-photon index diagrams, and study the 
distributions of the flare parameters. 

The sample is hosted by 5 blazars, of which the most ac- 
tive were 3C 454.3 (16 flares) and PKS 1510-089 (14 flares). 
In general, flares produced by 3C 454.3 are longer and more 
complex than those produced by other sources. This can 
be attributed to a very small viewing angle of 3C 454.3 
(dobs "C #j), a consequence of which is that a larger frac- 
tion of the jet volume is strongly Doppler-boosted. Other 
blazars, having o ^ B ~ 0j, produce shorter flares that are 
time-asymmetric, with the flare onset shorter than the flare 
decay. This is because jet regions that are more Doppler- 
boosted are located closer to the observer, and thus brighter 
subflares precede fainter ones. Long multiply-peaked flares 
indicate fragmentation of the gamma-ray emitting regions 
in blazar jets. I am unable to identify any regularities in the 
spectral variations of flaring blazars. 
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